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^ Abstract. We calculate the conditions required to produce a large local trispectrum during 

c/5 two-field slow-roll inflation. This is done by extending and simplifying the 'heat-map' ap- 

, ^, proach developed by Byrnes et al. The conditions required to generate a large trispectrum are 

broadly the same as those that can produce a large bispectrum. We derive a simple relation 

►^ between tnl and /nl for models with separable potentials, and furthermore show that ^nl 

■^ and Tnl can be related in specific circumstances. Additionally, we interpret the heatmaps 

^sl" dynamically, showing how they can be used as qualitative tools to understand the evolution 

of non-Gaussianity during inflation. We also show how /nl, t^nl and ^nl are sourced by 

^^_^ generic shapes in the inflationary potential, namely ridges, valleys and inflection points. 
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1 Introduction 

Observational constraints on the statistics of the primordial curvature perturbation provide 
a powerful test of inflationary models. For models driven by canonical scalar fields, the 
distribution of the curvature perturbation is extremely close to Gaussian at horizon crossing 
[1-4]. Moreover, after horizon crossing, the curvature perturbation is conserved for single- 
field models [5-11]. In multi-field models, however, if isocurvature modes play a role, it 
can evolve, and this evolution must be followed (for example see [12-22]). As a result the 
statistics may become less Gaussian, and even sufficiently non-Gaussian to be detectable by 
future probes, such as the Planck satellite [23]. In this paper, we consider models in which the 
non-Gaussianity of the curvature perturbation can evolve to become large during a slow-roll 
inflationary phase [24-27]. Moreover, we focus on two-field separable potentials, for which 
analytic formulae are available for the non-Gaussianity parameters [25, 28-31]. 

The best constraints on local non-Gaussianity to date are provided by the analysis 
of data from the WMAP satellite. The bispectrum is parametrised by /nLj which is con- 
strained as —5 < /nl < 59 at 95% C.L. [32]. tnl and ^fNL parametrise the trispectrum and 
are constrained as -0.6 < tnl/10'' < 3.3 and -7.4 < g^h/W^ < 8.2 [33, 34], with Ref. [35] 
finding the slightly different constraint —5.4 < 5nl/10^ < 8.6. Planck will improve on these 
constraints considerably, and in the absence of a detection is expected to give the bounds 
|/nl| < 5, Tnl < 560 [36] and |(7nl| < 1-6 x 10^ [33]. If a given inflationary model predicts a 
magnitude of the /nl, ''"nl or gjstL parameters greater than these forecast bounds, then we re- 
fer to such a prediction as 'observationally relevant'. This is in contrast to a non-Gaussianity 
that could in principle be measured by an ideal observation, which conservatively requires 
these parameters to have magnitudes greater than unity. Following the standard diction in 
the literature, we describe such models as producing a 'large' non-Gaussianity. 

Our aims in this paper are as follows: First, to show the constraints on the initial 
and final conditions of the inflationary evolution that lead to large values of the (7nl and 
Tnl parameters, and to compare these to the constraints that are required to produce a 
large value of /nl- Secondly, to derive a simple relation between tnl and /nl for two- 
field slow- roll models with separable potentials, which is compatible with the well-known 
consistency relation [37]. We will also investigate under which conditions it is possible to 
deduce robust relations between (7nl and tnl- Thirdly, following earlier studies of the /nl 
parameter [38, 39], we discuss the behaviour of tnl and gfNL associated with generic features in 
infiationary potentials, namely ridges, valleys and infiection points. Finally, we demonstrate 
that, in addition to showing the spectrum of possible behaviours, our formalism can be used 
as a qualitative tool to give insight into how the non-Gaussianity parameters evolve during 
inflation. 

Our method is an extension of the graphical approach employed by Byrnes et al. [27] to 
study the /nl parameter. The method begins by re-casting the expressions for /nl, tnl and 
5nl as a sum of terms. Each term is the product of a combination of slow-roll parameters and 
a function of two other parameters related to the initial and flnal conditions satisfled by a 
given evolution. Since the combination of slow-roll parameters is guaranteed to be much less 
than unity, if the non-Gaussianity parameters are to be large, the functions must take values 
much greater than unity. One can therefore make 2D contour plots (or 'heatmaps') of these 
functions, and so identify the regions in parameter space which can give rise to large local 
non-Gaussianities. In this study, we have found that it is possible to simplify the expressions 



to such an extent that only one heat-map is required to understand the conditions for a large 
bi-spectrum, one additional heat-map is required for tnl, and three further heatmaps are 
required for ^nl- 

Recently, Peterson & Tegmark have also undertaken a study of the bi- and tri-spectrum 
parameters [40-42] in the setting of slow-roll inflation, arriving at compact relationships 
between tnl, 5nl, /nl and the tilts of the curvature and isocurvature power spectra. Our 
approach is complimentary to their study, and provides new insights. For example, working 
directly with analytic expressions for /nl, tnl and ^nl allows us to present more explicit 
formulae, and in particular helps us to understand the evolution of non-Gaussianity directly 
in terms of the shape of the inflationary potential. Our general conclusions are, however, 
very similar. 

The paper is structured as follows: §2 provides some background material and reviews 
how analytic expressions for observable parameters can be derived using the 5N formalism. 
The expressions themselves are given in appendix A. §3 shows how we can simplify these 
expressions for the bi- and tri-spectra which we then analyse in §4.1 and §4.2 respectively. 
This analysis motivates three generic shapes in the inflationary potential which we study in 
§5, finding the peak values of /nl, tnl and (7nl in each case. In §6 we give concrete examples 
which illustrate the findings of this paper. We conclude in §7. 

2 Background Theory 

We consider inflation driven by two canonical scalar fields (pi with i = 1,2, self-interacting 
through a potential W{4'i, 4>2)- The scalar field equations of motion are 

4>i + 3H^i + W,i = 0, (2.1) 

where a comma denotes partial differentiation with respect to the field (pi , dots are derivatives 
with respect to proper time t, and H is the Hubble rate given by the Friedman equation 

i 

We choose definitions of the slow-roll parameters as 
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tt (2.3) 



rjij = M^,-^, il, = Ml,V2e 



Expressions for ^?-^ and rjij in the frame adapted to the field evolution, which we refer to as 
the kinematic frame, are given in appendix A. Inflation occurs when e < 1, and the so called 
'slow-roll limit' is given by e <C 1. In this limit the fields' kinetic energy can be neglected, 
and the field equations are well approximated by the slow-roll equations of motion 

2,H^i + W^i = {), SM^iH"^ = W . (2.4) 



2.1 The Curvature Perturbation 

The primordial curvature perturbation on uniform density spatial hypersurfaces (see for 
example [43]) is denoted by (. 

The statistical properties of C, which are constrained by observations, are commonly 
measured in terms of its power spectrum, bispectrum, and trispectrum, defined as 

(Cki Ck2 Cka) = {2TTf5\ki + k2 + k3)Sc(^l' ^2, A:3) , 
(Cki Ck2 Ck3 Ck4) = (27r)353(ki + ka) + k3 + k4)rc(fci, A:2, h^ h) , (2.5) 

respectively. For the local shape of non-Gaussianity, it is convenient to further parametrise 
the bispectrum and trispectrum in terms of the dimensionless parameters, /nl, ''"nl and gjstL, 
defined by 

Bcih, k2, h) = -hL [Pdki)Pc{k2) + 3 perms] , (2.6) 

T(;iki,k2,k3,k4) = TNL [P(;{ki3)P(^{k3)P{k4) + 11 perms] 

54 

+^5NL [Pdk2)Pc{k2)Pdh) + 3 perms] . (2.7) 

2.2 The 6N formaUsm 

In order to follow the evolution of C on super-horizon scales, and calculate its statistics, we 
employ the separate universe approach to perturbation theory [6, 44], and the 6N formalism 
[14, 28, 45]. In this approach, spatial gradients are neglected on scales greater than the horizon 
size, and each spatial point is assumed to evolve as a separate FRW universe. Choosing a 
flat initial slicing at the time at which observational scales crossed the cosmological horizon, 
t = t* , and a later uniform density (constant H) slicing at t = to then ( on the final slicing 
can be equated with the excess expansion, C = 6N. 

During slow-roll inflation, field velocities are functions of field positions. Taking this 
to be a good approximation at horizon crossing, the subsequent number of efolds undergone 
by any 'separate universe' is a function purely of the initial field values on the flat slicing, 
N{(j)l, (/jg), even if the evolution subsequently evolves away from slow-roll. One therefore finds 
that 

C = 5N = NM\ + j^N,,5cl^*M* + j^Ni,k5(t>*5<i)*54>l + . . . , (2.8) 

where here and throughout we employ the summation convention, N is the number of efolds 
from t* to tc, a subscript i on N represents a derivative with respect to the light fields 
at horizon crossing (/>^, and 54)* are the field fluctuations on a flat hypersurface at horizon 
crossing. 

Using the expression for (^ from Eq. (2.8), and recalling Eqs. (2.5)-(2.7), one can write 
observational parameters in terms of the derivatives of N [14, 28-30], 

f - iM^N^f, , n. 1 - ^^ (iVfciVfc)2 ^' ' ^^-^^ 
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where n^ is the spectral index, P* = //^/(27r)^, and where we have assumed the field fluctu- 
ations to be Gaussian at horizon crossing, which we recall is an excellent approximation for 
canonical fields [2-4]. 

2.3 Analytic formulae 

To calculate the observational parameters for a given inflation model we require the deriva- 
tives of N. These can always be calculated numerically, but analytic progress is possible only 
when the slow-roll equations of motion, Eq. (2.4), are a good approximation, and when a spe- 
cial 'separable' form of potential is assumed. This can either be a sum-separable potential of 
the form VF = Ei Vi{4>i) [15, 25, 29, 46], or a product-separable form, W = IliVi{(t)i) [47]. It is 
also possible to extend the analytic analysis slightly to models of the form W = (^^ Vi(0j))^' 
[48], where A is an arbitrary constant. For simplicity, we will restrict our study to two-field 
sum- and product-separable potentials, labelling the fields (pi = (j) and 02 = X- The two-field 
analytic formulae we employ for the derivatives of A^, and for the observation parameters 
which are derived from them, are given in appendix A. 

2.4 The evolution of statistics and the adiabatic limit 

During inflation, if a given field space evolution follows a straight trajectory, then the statistics 
of the curvature perturbation will remain constant. Conversely, if the trajectory turns, the 
statistics will evolve [18, 24, 25]. During a turn, typically in the early stages, it is possible 
for the statistics of the curvature perturbation to become highly non-Gaussian [26, 27, 49]. 
If inflation were then to end suddenly, perhaps through a hybrid transition, it is possible 
that this large non-Gaussianity is preserved into the subsequent phases of the universe's 
evolution. Often, however, if the turn is completed, the non-Gaussianity returns to negligible 
levels. This possibility was emphasised in Ref. [31, 50]. 

This is not the only possibility, however. It is also possible for the bispectrum [38, 39, 
51, 52] and trispectrum [53] to be large once a turn is completed, even after all isocurvature 
perturbations have decayed. If isocurvature modes decay, ( and its statistics are subsequently 
conserved on super-horizon scales [7, 8], and this has recently been referred to as the adiabatic 
limit [39, 54]. As discussed extensively in the literature (for example see [13, 15, 17, 39]), 
reaching such a limit greatly simplifles the task of making observational predictions for an 
inflationary model. We emphasis, however, that there is no requirement that this occur, 
and if it does not, isocurvature modes will persist at the end of inflation [55]. The decay of 
isocurvature modes during slow-roll inflation occurs if the evolution reaches a regime which 
is effectively single field, such as a valley region with steep sides which force isocurvature 
perturbations to decay. One purpose of our study is to systematically identify all of the 
conditions which lead to a large tnl and ^nl during slow-roll infiation, both in cases where 
isocurvature modes have, or have not, decayed. 

If the evolution reaches a straight line trajectory along one of the axes, or the isocurva- 
ture modes decay, there is a considerable simplification of the formulae for /nl, 5nl and tnl 
given in appendix A. In particular, the terms involving Z, A and B^ (or Ap and Bp ) as de- 
fined in Eqs. (A. 2), (A. 6) and (A. 12) (or eqs. (A. 21) and (A. 27)) tend to zero"*^. These terms 

^As discussed at length elsewhere [38], there is a possibility that, for sum-separable potentials, Z does not 
go to zero. This exception only occurs if one of the fields, (f)k, is completely orthogonal to the final straight 
line direction, and Vfc((/>fc) tends to a constant as Vkitpk)' tends to zero. Such a possibility can be avoided 
since we are free to reparametrise the potential and associate the problematic constant with another of the 
fields. 



vanishing is referred to as the Horizon Crossing Approximation (HCA) [56]. This greatly 
simphfies the expressions (A. 14) and (A. 28), aUowing us to make much stronger statements 
about the relations between the observable non-Gaussianity parameters. 

A final clarification is in order. Throughout this paper we calculate statistics on uniform 
density hypersurfaces. If, however, inflation ends suddenly, a further complication is that the 
surface on which inflation ends may not be a uniform density one. If this is not the case, 
an additional source of ( will be produced [57], which will alter the statistics we calculate. 
We do not discuss this possibility further here, except briefly in §6.1. We note that in cases 
where all isocurvature modes decay there is no such possible additional effect. 

3 Analytic non-Gaussianity 

In this section we present the central results of our paper. We reformulate the analytic 
formulae for the non-Gaussianity parameters, which are reviewed in appendix A, into a set 
of terms, each of which is the product of a slow-roll suppressed part, and a function of two 
variables related to the initial and final conditions. Only the functions can be large, and so by 
plotting them as two-dimensional heatmaps, we can identify the initial and final conditions 
which lead to large non-Gaussianities. A key part of this procedure involves neglecting 
terms in these expressions that are always too small to contribute significantly to a large 
non-Gaussianity. We will see in subsequent sections, that these heatmaps are also useful for 
understanding the dynamical evolution of the non-Gaussianity parameters. Byrnes et al. [27] 
used a similar method to study the /nl parameter. We first improve this bispectrum analysis 
by reducing to one the number of relevant heatmaps, and then consider the trispectrum. 

3.1 Variables 

We begin by defining the first important variable, the angle 9 in the {(j), x} phase space 
in which the infiationary trajectory evolves. This can be done in terms of the slow- roll 
parameters e^ and e^^, defined in Eq. (2.3), as 

!^ = cos2 0, ^ = sin2 0. (3.1) 

Since we assume that both fields are monotonically decreasing (which follows from our use 
of the slow-roll equations of motion), 6 is constrained to lie in the range < ^ < 7r/2. 

The 5N expressions, given in appendix A, also involve the quantities u and v, defined 
in Eqs. (A.l) and (A. 17), which lie in the range zero to one. It proves convenient, therefore, 
to define a second angle, a, in terms of these quantities as 

u = cos^ a , V = sin^ a. (3-2) 

We note that in the product-separable case, a = 9. The situation is not so simple in the 
sum-separable case, as we shall discuss in §4.1. 

Substituting these definitions into the expressions for the non-Gaussianity parameters, 
Eqs. (A. 16) and (A. 30), we can eliminate the variables u,v,e^ and e^ in favour of a, 9 and 
e. The observables are then only functions of a, 9 and 9* , as well as e and the other slow-roll 
parameters. Some of these slow-roll parameters are evaluated at horizon crossing, whilst 
others are evaluated on a later constant energy density hypersurface, usually labelled 'c'. In 
the following we can drop this label without ambiguity, since within the 6N formalism all 
quantities evaluated after horizon crossing are calculated on a uniform density hypersurface. 



Quantities without a '*' label are therefore assumed to be calculated on this later-time 
uniform density hypersurface. 

3.2 Bispectrum 

For the /nl parameter the procedure outlined above leads to the expressions 
6 



r, /nl = /le* - f2V*ss + hvls + 2/ ri {riss - e) - Sum separable , 

5 

^/nl = -firfss + fsvls + 2/??ss - Product separable , 

5 



(3.3) 
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where, similarly to Ref. [27], we have defined the functions f{a,9*), fi{a,9*), A and (7 as 

/ = ^^2 (^°^ a-cos 6' ) , /2 

^ sin^ 29* 

A = cos a sin 9* + sin a cos 9* , Q. — no 

W^ sin"' 2a 

(3.4) 

Since < il < 1 ^, the functions /i-s-s and / all multiply quantities of 0{e) or smaller, 

and so a necessary, though not sufficient, condition for /nl to be large is that the magnitude of 

one or more of these functions is large. Our analysis, therefore, identifies only the conditions 

for which it is possible to produce a large /nl during slow-roll inflation. 

Because we are only interested in cases where /nl can be large, the expressions (3.3) 
may be further simplified by noting that | /i | is bounded by order of unity and so the term 
fi€* is negligible. Similarly, the term fsr]*^ is negligible and can be dropped. This latter 
result follows by recalling that we have the freedom to interchange between rj^u, rj^s and rjss 
via the relations 



r](js = 5 tan 29 {rjss — rjacr) - Sum separable , 

rjas = 5 tan 29 {rjss — ijacr + 2e) - Product separable . 

Considering first the sum-separable case and expanding /s?]*^ as 



(3.6) 



/3<,= [sin2r/3]7?:,+ [2(l-sin2r)tan2r/3](r?:,-r/:j, (3.7) 

our result follows from the fact the the terms in square brackets can never become larger than 
order of unity. This procedure works because rj*^ tends to zero in the limits of 9* — )• 0, 7r/2 
which counters the divergence in f^. Since Eq. (3.6) has the same form for both sum and 
product-separable potentials, we may neglect the term /st/*^ in both cases. 



2 



This follows by using (A.l) and the associated definitions found in appendix A to find 

[U-e + Ve^ - [/6x)(V*e - Ve^ + Ue^) " {Ue + Ve^ - Ut-^){Ve - Vt^ + Ue^) ' ^ ' ' 

where the second inequality follows using U < U* and V < V* . 



A final simplification follows by noting that /2 = 1 + /• We arrive at the extremely 
simple approximate expressions for /nl 
6 
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Vls + '^^iVs 



:/nL ^ f -V*ss + 2f?s 



- Sum separable , 

(3.8) 

- Product separable , 



5- 

where / is positive definite. These simpler expressions make it transparent that the condition 
for /nl to be large is that / » 1. We emphasise that these approximate expressions will be 
extremely accurate when /nl > 1- 

3.3 Trispectrum 

A similar analysis can be performed for the tnl and ^nl parameters. For sum-separable 
potentials we find 

tnl = nr]*/ - T2r]*^r]*^ + t^t]*^/ - ue*r]*^ + rge*??;^ + Tee*^ 

-T7 ^ r]*^{r]ss -e) - ts ^i ??* ^ (?7ss - e) + TgQe*{r]ss - e) + 4:Tn^ {rjss - ef , (3.9) 

54 o o 1 1 lolo 

3 3 

+91 ^""^ {ilss - '^VasiVss + e)) +Ag3n — cos 29 rj^sivss - e) , (3.10) 

and for the product-separable potentials we find 

TNL = nry*^^ - T2r]l,rils + ^3^/*,^ - T^V|l^^^ss - rsvlsVss + ^ttjI^ , (3.11) 

3 3 

-^'^V*ssVss - -Tsril.riss + gi {fsss - '^VcTsVss) + 4fif2?7L > (3-12) 

where the various functions occurring in these expressions are defined as 

n = -^(cos^ a sin^ 6* + sin^ a cos^ 0*) , T4 = 2/1/2 , 

A'^ 

r2 = — 7^5 — (cos asm t) — sm a cos t' ) , T5 = z/1/3 , 



(3.13) 



/l / 8 ■ 2 /)* , ■ 8 2 /)*\ ^2 

T3 = — "2 (cos a sm + sm a cos U ) , tq = Ji , 

Ty = — (cos a sin 9* + sin a cos 6**) , rg = 4/i/, 

sin 20* sin^ 2a / 9 2 /i*\ /sin^2a 

^8 = ^ (cos a - cos ) , T = ^^ , 

gi = g3sm2a, 52 = 53 cos 2a, 53 = - — (cos^ a - cos^ 61* ) . 

We can now proceed similarly as we did for the bispectrum and show that a number 
of these terms are negligible if the trispectrum parameters are large. Since the trispectrum 
functions pre-multiply quantities that are second order in slow-roll, in this case we neglect 
any functions that are never larger than 10, rather than order of unity. The details of this 
analysis can be found in appendix B. 
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(3.14) 



(3.15) 

A further simplification can be made by noting that C > 1 and that /i is at most of order 
unity. The second terms will therefore be suppressed relative to the first by at least 0{e*). 
We thus find 

6 -' 



TNL^^CI -/nl) , (3.16) 

which is valid for both sum and product-separable potentials. 

In ref. [33], the ratio between tnl and (|/nl) was parameterised by Anl- Peterson et 
al. subsequently showed that ^nl = l/'^c foi' two-field models under slow-roll [41], where re 
determines the fraction of the curvature perturbation which is sourced by a horizon crossing 
isocurvature mode. Our result, Eq. (3.16), is complementary to this analysis, explicitly 
showing the form of r^. in terms of the dynamics of infiation for separable potentials. 

3.3.2 5NL 

Unsurprisingly, g'NL does not simplify so neatly. In the product-separable case we find 

^5NL ^ TNL ( '" r, ^'^ ) - f/nl(2?/L + r]ss) - OiCsss"^ + 91 Csss " "^VcrsVss , (3.17) 
OU \Vss " ^"Hss J L -I 

where gi = \ (ts sin 2^* cos 29* — T2). For sum-separable potentials we obtain two additional 
terms, and the expression takes the form 

54 f ri*,-n {r]ss - e) \ 6 2 

— 5NL ^ TNL — ^TVT r - 7/NL(2r/^s + n {r]ss - e)) - Msss 



50"^^" ''^yr,ts-2n{7]ss-e)J 5 

+51 ^3/2 [eL. - 2r^UVss + e)] - lfife*r]*ss 

f W \ 

+Ags n {r]ss - e) ( :rr^ COS 2%s - ^ cos 2a{r]ss " «) j • (3-18) 

These formulae for q'nl are once again complimentary to those derived by Peterson 
et al. [41]. One advantage of our results is that because of their explicit nature, we can 
easily consider which shapes in the infiationary potential will cause the different terms in 
Eqs. (3.17) and (3.18) to become large. Furthermore, the Horizon Crossing Approximation 
is easily implemented in our formulae by taking [7—7-0. 

4 Analysis 

We now turn to the interpretation of the formulae we have presented in the previous section. 



4.1 Bispectrum 

The expressions (3.8) imply that a necessary condition for |/nl| > 1 is that / ^ 1. We can 
see when this occurs by plotting /(a, 9*), as shown in Fig. 1. 
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Figure 1. Hcatmap of the function /. Since /nli Vss a-nd e are symmetric under field exchange 
(/) -f-)- X, the function / must have the same symmetry. We can see this by inverting the heatmap 
though the point {0*,a) — (7r/4, 7r/4) which leaves / invariant. 

For product-separable potentials, a and 9 are trivially related as a = ^. The initial 
(horizon crossing) value of a for product-separable potentials is therefore ajnit = 9*. For 
sum-separable potentials the same initial condition is found. This can be seen by taking the 
limit c — 7- * in Eq. (A.f). This means that all separable inflationary models start on the 
diagonal of the heatmaps at horizon crossing. Since / = on the diagonal, we know that 
initially /nl will be given by the various negligable terms that we have dropped in (3.8), and 
so |/nl| ^ 1 initially, as we expect [2]. 

As the model evolves from a given 9* , a varies such that the model traces a vertical line 
on the heatmap. For this trajectory to ever intercept one of the regions in which / 3> f the 
initial conditions must be such that the initial phase-space velocity is dominated by one of 
the two fields (/> or x- 

4.1.1 The product-separable case 

As an example, let us take a product-separable model with 9* <ti 1, such that the horizon 
crossing conditions correspond to a position on the diagonal of the heatmap of / in the lower 
left-hand corner. If the angle 9 (= a) increases towards tt/2, there can be a transient 'spike' 
in /nl as the model passes through the region in which / is large. If the trajectory turns back, 
so that 9 decreases, it may well pass back through this region again and another transient 
signal in /nl can be produced. Of course, whether or not a significant spike will occur is also 
dependent on the slow-roll parameters that define the model. The magnitude of this spike 
in /nl increases as 9* decreases towards zero. We must be careful with the interpretation, 
however, since our analysis presumes separable potentials and this places strong constraints 
on the possible modes of behaviour. If the initial field velocity is exactly aligned to one of 
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the field axes then it wih (classically) remain so indefinitely, since this implies 6 = 6* always. 
Furthermore, for neighbouring initial conditions with sufficiently small values of 6* , one may 
find that 6 will take longer than the required length of inflation (roughly 60 efolds) to grow 
sufficiently for there to be a significant enhancement of /nl • There will therefore be an upper 
bound on the value /nl achievable by any such potential'^. 

4.1.2 The sum-separable case 

To make progress with understanding the dynamics in the sum-separable case, it is necessary 
to understand the relationship between a and 6. To do this we differentiate Eq. (3.2) to find 

W*sin^2a ^ ' 

where a prime represents a derivative with respect to some time variable. The fractions in 
Eq. (4.1) are positive definite and so a increases as 6 increases and vice versa. 

Vacuum domination. More progress is possible by considering the vacuum-dominated 
limit, where W ~ W* . This is a good approximation for some models of inflation such as 
hybrid inflation. Eliminating the ratio W/W* from Eq. (4.1), we see that a c:^ 6. Moreover, 
one also finds fi ~ 1 in this limit. Consequently, for vacuum-dominated sum-separable 
potentials, one may use the heatmaps in the same way as for product-separable potentials. 
Furthermore, we see that in this case the sum and product-separable formulae for /nl are 
identical, apart from the presence of the slow-roll parameter e in the sum-separable case. 
Thus for models with r]ss ^ e there is a very near equivalence between product-separable 
potentials and vacuum- dominated sum-separable potentials in terms of their contribution to 
/nl- 

General models. We now consider Eq. (4.1) for general sum-separable models without 
vacuum-domination. The ratio W/W* is initially unity and decreases towards zero. This 
ensures that, whilst the angle is monotonically varying, a lags behind 6. Furthermore, the 
difference between a and 6 will become more pronounced the smaller the ratio W/W* . When 
this ratio goes to zero, we see that a remains constant despite any subsequent turning of the 
trajectory in phase space. Similarly, a will cease to evolve if either of the limits — )• or 
6 — > 7r/2 are reached. Physically, these limits correspond to straight lines in the field phase 
space under which conditions it is well known that (" does not evolve. We note that Q is zero 
in any of these three limits and so we are only left with the 7]*^ term in /nl- This is the 
approximate HCA formula for /nl- We note that since a j^ 6, it is quite possible for / to be 
large when a becomes a constant, and so produce an observationally relevant constant /nl- 
For a given evolution, the final value of a is readily extracted once the initial and final field 
values are known, and hence one can check whether the correct value of a is reached in order 
to give a large non-Gaussianity. 

The final case to consider is when the trajectory turns back on itself during its evolution. 
There is no barrier to constructing sum-separable potentials which exhibit this behaviour for 
particular evolutions. This will mean that the model moves up and then down a vertical line 
on the heatmap, perhaps many times. To fully quantitatively understand how this movement 

^There will also be another upper bound on J/nlI imposed by the quantum scatter of the field near horizon 
crossing which will prohibit the limit S* — > being physically realisable. 
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occurs requires a knowledge of how the potential W varies with 9, which is necessarily model- 
specific. However, it is possible to gain some more detailed intuition for how a varies with 9 
by rewriting Eq. (4.1) as 

h'{a) = ^h'{9), (4.2) 

where we have defined the function h(x) = 4a; — sin 4^;, which is monotonic in x. One sees that 
h{a) and h{9) increase and decrease simultaneously and hence the velocity h'{a) is always 
smaller than the velocity h'{9). Since W/W* is constantly decreasing then so is the range of 
values of h{a) which the evolution can reach. Ultimately h'{a) — ?• and h{a) takes a constant 
value. Since h{a) is monotonic in a we see that restricting the range of h{a) translates into 
introducing 'excluded regions' at the top and bottom of the heatmaps, the size of which will 
grow as the potential drops, and ultimately the whole of the heatmap will be excluded except 
for the final value of a. Once again, for a given evolution the final value of a can be readily 
calculated, and one can check if it is in the observationally relevant regime. 

4.1.3 The role of n 

The value of fi determines which of the terms r]*^ or r]ss dominates in Eq. (3.8). For product- 
separable or vacuum-dominated sum-separable potentials, ri is effectively fixed to be unity, 
and so which term dominates depends on how r/^s evolves during the evolution. On the other 
hand, if il — )• then the r]*g term will dominate. It is instructive, therefore, to think of 
inflationary models belonging to one of two classes, those for which il ~ 1 throughout and 
those for which fi — )■ at some point. We now briefly consider each of these cases in turn. 

Evolutions for ■which i7 ~ 1. These models either have product-separable or vacuum- 
dominated sum-separable potentials. In the simplest cases of interest, such as falling from 
a potential ridge or rolling into a vacuum dominated valley [38], the absolute value of the 
potential does not change significantly during this phase of evolution. If rjss — vts then one 
finds |/nl — fvts- Since / is positive definite, we see that the sign of /nl is the same as 
the sign of the isocurvature mass. This is in agreement with previous calculations [38] where 
a ridge shape {r]ss < 0) was shown to lead to a negative /nl and a valley shape {r]ss > 0) 
was shown to lead to a positive /nl- If inflation does not end abruptly, for example by a 
hybrid transition, but the fields continue their evolution, then /nl will continue to evolve 
until — )• or — 7- tt/2. In either of these limits / — )• and so /nl is much smaller than 
unity. 

Evolutions for which ^2 — )• 0. These models naturally reach a limit in which a becomes 
constant, and the HCA becomes a good approximation. This may or may not coincide with 
a scenario in which the isocurvature perturbations decay and an adiabatic limit is reached. If 
an adiabatic limit is arrived at, then we can be certain that there will be no further evolution 
of C, and its statistics. When i7 — )• we find |/nl — —fvts- Due to the minus sign, if we 
begin in a region with a negative isocurvature mass, for example on a ridge, then such a 
model will eventually produce a positive /nLi whereas if we begin in region with a positive 
isocurvature mass, then /nl will ultimately reach a negative limiting value. 

In summary, the positive /nl that is marginally preferred by the WMAP data [32] can 
be generated in two possible ways. Either, the evolution must begin with a small 9* (or one 
close to it/2) and an initially negative isocurvature mass, and then evolve until a naturally 
takes a small constant value suitable to make / large. Or alternatively, the evolution must 
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evolve from a small 9* (or one close to vr/2) during which time the isocurvature mass is 
positive, and some mechanism much interrupt the dynamics whilst / is large. 

4.2 Trispectrum 

We now turn to analyse the trispectrum. Our aims are twofold. First, to understand the types 
of models and the initial and final conditions for which the observational non-Gaussianity 
parameters tnl and g'NL can be large. Secondly, to infer if it is possible to relate these non- 
Gaussianity parameters to one another, perhaps for specific classes of models. Such relations 
are potentially important in order to observationally exclude classes of models. 

Referring to Eqs. (3.14) we see that tnl can only be large when r or / are large. We plot 
these functions side-by-side in Fig. 2. The fact that both functions peak in similar regions of 
the parameter space means that the classes of models that are capable of producing a large 
Tnl are the same as the classes of models that can produce a large /nl- 
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Figure 2. Heatmaps of / and r on the same scale. We note that both are large in the same regions 
of the parameter space. 



4.2.1 Tnl and /nl relation 

For canonical scalar field inflation, tnl and /nl satisfy the condition [37] 

2 



"^NL > ( ^/nL 



(4.3) 



where equality occurs for single field infiation'"^. This is fully consistent with our result (3.16), 
once we recall that C > 1. We plot C in Fig. 3. The interesting regions of this plot are those 
for which tnl is not small and so /^ ^ 1 as shown in Fig. 2. We see that if a model has 
9* <^ 1 and 9 subsequently increases, then such a model will first enter the region in which 
C S> 1 and so tnl will grow whilst /nl remains small. It is unsurprising that tnl evolves 



*See also the recent work by Sugiyama [58] claiming that this equality is broken when calculations contri- 
butions are included from all loops. 
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first, since being associated with a higlier order moment, it will be more sensitive to outliers 
of the 6N distribution. For large a, one can see that C ~ 1 to a very good degree of accuracy, 
and the single field relation becomes a good approximation. We note that tnl deviates more 
from the single- field limit C = 1 when 6* is fine-tuned to be closer to the 6* = axis. 




Figure 3. Heatmap of the function C illustrating the conditions required for a model to deviate from 
the single field consistency result for tnl- This is only physically interesting in the regions for which / 
is large which are close to the sides of the heatmap. The region where C ^ 1 overlaps with part — but 
not all of — the region where / is large and one can see that this overlapping region is displaced from 
the sides of the heatmap. 



4.2.2 (7nl and tnl relations 

We now turn our attention to consider possible relations between /^nl and tnl- 

The rj = case. First, we consider the case of non-vacuum dominated sum-separable 
potentials for evolutions which reach f] = such that all of the observational parameters 
have ceased to evolve. From Eq. (3.18) we see that 



54 

50 



5NL ^ TNL - -t/nl (4r/*5 - he*) - g^C 



(4.4) 



The /nl term in Eq. (4.4) is suppressed relative to the tnl term by a relative factor of 
/nl X £)(€*) and so may be safely neglected. In cases were ^*^^ can be neglected, for example 
in the absence of significant terms in the potential beyond quadratic order at horizon crossing, 
we find a relation between tnl and ^nl as 



54 
50 



ffNL - TNL • 



(4.5) 



Next, let us look at the term g^Qss and assess when it may be relevant. If an inflationary 
potential exists for which this term is important then it will have a signature such that ^nl 
deviates from Eq. (4.5). The heatmap for (74 is plotted in Fig. 4 and we see that the areas in 
which 54 is large are very small in comparison to the corresponding areas for the r function. 
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Practically speaking, this ensures that one has to tune the parameters of the model to a very 
high degree in order to access this region. In §5, we show that, in the interesting limit where 
6* and a are small, we can accurately approximate 2g4 ~ 6*t. Thus a necessary (but not 
sufficient) condition for I^nlI to be large is that 9*t S> 1. In addition, for g^i^ to deviate from 
Eq. (4.5) we require that L = 9*^*^^ /'^ffss is not small. Potentials with L ^ 1 therefore 
have the capacity to generate |(7nl| ^ "^nl- 

Considering a sum separable potential, in the case where 0* ^ 1, to leading order we 
find 



% 



* 2 



£* 






4^1 



(4.6) 



where we have indicated r\^^ = r]^ and £^^^ = £^. Expanding out the slow- roll parameters 
in terms of potential derivatives, one finds that L ~ ^,x^,xxx/^^,yx' ^^ °^^ considers a 
general power-law potential V{x) = Vqx" then L = for n = 0, 1, 2 and L < 1/2 for n > 3. 
On the other hand, for an exponential potential V{x) = Vqc^^ then L = 1/2. In these two 
cases we would therefore not expect to find a deviation from Eq. (4.5) beyond a factor of 2. 
It is interesting to note that a potential of the form V{x) = a ln(x — b) for constants a and b 
has L = 1 and so for such a potential the two leading order terms in (^nl exactly cancel and 
so 5NL is of order /nl xO{e*). 

We now ask if there is a potential with L ^ 1. Considering a polynomial potential V{x)i 
the necessary condition is for the potential to possess a linear term, a negligible quadratic 
term and at least one term beyond quadratic order. The simplest such potential is a sloping 
inflection point of the form V{x) = Vq + hx + q^X^- Nearby the inflection point one has 
r/jj, ~ whilst ^ Cv = ^pi^-^/^o ^^'^ ^^ ^ diverges. We shall examine such an inflection point 
further in 56. 
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Figure 4. Heatniap of the function g^ which is appropriately antisymmetric about field exchange. 
The very small lobe in the bottom left hand corner has positive 54, whilst the lobe in the top right 
hand corner has negative 34. This pattern is repeated for gi and g^ as seen in Fig. 5. 

The ri 7^ case: It is considerably harder to make concrete statements about the value 
of the trispectrum parameter /^nl, and its relation to tnli when Q 7^ 0. We can take a step 
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in this direction by plotting the heatniaps for the remaining functions which appear in the 
expressions (3.17) and (3.18), gi and §3. These are shown in Fig. 5. The similarity between 
these plots, and those of / and r, tells us that the types of inflationary potential, and initial 
conditions, that can give rise to a large g'NL a^^e similar to those that may give rise to a large 
/nl and tnl- No new regions of interest appear for (^nl which were not present before. 
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Figure 5. Plots of gi and (73 on the same scale as used for the other heatmaps. We note that they 
are very similar and are anti-symmetric under field interchange. Both are large in the same regions 
of the parameter space as t. 

The ri ~ 1 case: Finally, we consider the case $7 ~ 1. This arises for the vacuum domi- 
nated potentials, and so can be relevant for models where inflation is terminated suddenly, 
perhaps through a waterfall transition. During a vacuum dominated phase of evolution, the 
flelds evolve only very slowly, and likewise turns in field space progress slowly. 

It is instructive to consider how the non-Gaussianity parameters given by Eqs. (3.8), 
(3.16) and (3.18) inay be approximated for $7 ~ 1, ^ 1 and a power-law potential of the 
form V{x) oc x" where we presume n to be a positive integer greater than or equal to 2. In 
the regime of interest where 9^6* (such that the functions / and r may be large) we require 
\x\ ^ \x*\- If ?T- = 2, such that the potential has a parabolic shape, then r]ss is approximately 
constant and from Eq. (3.18) for (^nl we see that the first term will be negligible in this 
case. Furthermore, if the potential is well-described by an expansion to quadratic order we 
see that the other terms in ^nl are also suppressed with respect to /nl and so g'NL is not 
large. This was the case found by Byrnes et al. [49] in their study of vacuum dominated 
quadratic potentials. 

However, for n > 3 the situation is different. The isocurvature mass r]ss is now a function 
of X and \r]ss\ will grow as |x| grows. In the regime of interest where |x| ^ |x*l we find that 
I^^NL — 5TNL from the first term in Eq. (3.18) and so it is clear that (/nl is not generally 
a vanishing quantity when the potential is described by terms beyond quadratic order. We 
note that the remaining terms in Eq. (3.18) for (7nl will mostly be small in this limit, but 
we will anticipate that the gi^sss term will contribute notably to the value of g'NL during this 
transient evolution. We shall look at the specific case of n = 3 in §5.3. 
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5 Typical shapes in the inflationary potential 

In this section we briefly show how the formulae that we have found for a large /nL) '^nl 
and gjsiL can be used to make quantitative predictions in certain generalised settings. We do 
this by studying how the bi- and tri-spectra transiently evolve as the field space trajectory 
passes over particular features in the inflationary potential. Motivated by the analyses in the 
previous sections, as well as other work [26, 39, 49], we consider potentials with one of the 
following three general features: 

• a ridge, 

• a valley, 

• an inflection point. 

Our analysis leads to simple expressions which give the sign and the peak values of the 
observable parameters generated by these features. By Taylor expanding a real potential 
about the initial conditions corresponding to a given inflationary potential, it is often possible 
to approximate the actual potential by one of these features over a range of field values. This 
can allow quantitative information about the predictions of models of inflation to be obtained 
readily without the need for detailed calculations [38]. 

The analyses of heatmaps show that large non-Gaussianity parameters require small 
angles in the parameter space which correspond to potentials for which the motion is highly 
aligned to one of the field directions, which we take to be (f) without loss of generality. We 
can then expand the functions such as / in the limit where 9* and a are small to find 

(^^ + ^2)2' ^-(04 + ^2)3' w-^; 

and the other parameters are related as gi ~ Or, 2g^ ~ r and 2g4 ~ 9*t. We now consider 
each feature in turn. 

5.1 Potentials v^ith a ridge 

The simplest possible ridge we can consider takes the form 

W = Wo + g(l)-^mlx\ (5.2) 

where g and m^ are taken to be positive and the Wq term dominates. Since the potential 
is vacuum dominated and fi ~ 1, the product-separable and sum-separable formulae for 
/nl, ''"nl and ^nl are identical, because vacuum-domination also implies that e <^ rjss- To 
be consistent with our sign convention we stipulate that x < such that W^i > and 
< 9 < 7r/2. The initial conditions are fine-tuned such that the field initially moves almost 
parallel to the top of the ridge, with ^ ^ 0. The isocurvature direction is therefore almost 
precisely the x direction, and rj^s ~ —m^^/Wo is a constant. This then gives 



6^ ^ * *2 54 U 



:/nl ^ fvl , T^h ^ rrf , ■^fi'NL ^ ^/nl ??t • (5.3) 



In this case ^nl is subdominant, even to /nlj and so we do not consider this further. The 
peak value of /nl and tnl can then be found by maximising the functions / and t as 9 varies. 
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They peak when 9'^ = a/30* and 9'^ = y/29* respectively and so 

cr /NLlmax — 16 0*' ^NLlmax — 27 0^ ' V"-*-^/ 

The sign of /nl is negative due to the sign of r/* , and the amphtude of /nl reaches its peak 
after tnl- The peak of tnl is only slightly larger than the square of | /NLlmax' but due to 
the difference in peaking times it is quite possible to have tnl S> | /nl I ■ 

5.2 Potentials vi^ith a valley 

The simplest possible vacuum-dominated valley takes the form 

W = Wo + ^m2</)2 + ^mlx^ , (5.5) 

where m^ and m^ are taken to be positive and the Wq term dominates. We again consider 
9* <^ 1, and assume that m^ ^ m^, so that the initial motion is almost exactly parallel to 
the (j) direction. Once again, therefore, the isocurvature direction is well approximated by 
the X direction, and r]ss = m^/Wo is again constant. This then gives expressions for /nl, 
Tnl and (Jnl which are identical to those in Eq. (5.3), but we note that the signs of /nl and 
5nl are now reversed. 

5.3 Potentials ^vith an inflection point 

The analysis of §4.2.2 with il ~ 1 illustrated how (/nl is negligible for inflationary potentials 
that are approximately quadratic, but that (7nl ~ ''"nl for potentials with cubic or higher 
order shapes. We now consider the simplest such shape, with an inflection point of the form 

W = Wo+g4>+\\x'- (5.6) 

o 

The regime of interest is close to but below the inflection point with |x| ^ 1 such that <C 1. 
The non-Gaussianity parameters are only large in the regimes where / or r are large, which 
requires \x\ 3> \x*\- In this limit we find that the leading order non-Gaussianity is given as 

6 54 5 

^/nl ^ 2/r7ss , TNL ^ 4rr/f^ , ^^nl ^ g^NL • (5.7) 

For a general power-law potential, the numerical relation between (7nl and tnl will take a 
different value. We note that ^nl simply follows the evolution of tnl and so we need only 
calculate the peak values of /nl and tnl- We find 

\ /nl|_ - -1-08^ ^ -nl|_ - l-48f|^ . (5.8) 

where these peaks respectively occur when 9"^ = a/S^* and 9^ = ^Jl?>/?)9* . 

6 Concrete Models 

In the preceding discussion we developed a number of tools for analysing a given two-field 
inflationary model. In this section we demonstrate the usefulness of these tools, by showing 
concrete examples of simple models which can produce large values of tnl and/or (/nl • 
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6.1 Two-field hybrid inflation 

The original hybrid model was introduced by Linde [59] as an effectively single field model 
of inflation, where the job of ending inflation was assigned to a second field, known as 
the waterfall field. The interesting feature of this set up is that inflation does not end 
gracefully, but rather through a sudden transition, which allows more freedom in constructing 
the potential that drives inflation. This model has been extended to the multi-field case, 
where two or more scalar fields drives inflation before the waterfall field ends it. In this 
section we consider a two-field hybrid model of infiation with the potential [28] 

W = Wo(l + ^V<t>^^ + ^r?xx') - (/<A' + 9<PX + hx^) y + ^^^ V;^ , (6.1) 

where the waterfall field ^ is held at during inflation, rolling rapidly to it's true minimum 
when 

ml = fct)^ + g4>X + hx^ • (6.2) 

Since inflation ends suddenly in this model before isocurvature modes decay, we need to 
account for the possibility that inflation might end on a surface other than one of constant 
energy density. This implies that the curvature perturbation produced by this model can be 
split into two parts, 

C = Cinf + Ce , (6.3) 

where Cinf is the part generated during inflation, and Ce is the part generated from the end 
of inflation dynamics. The condition (6.2) for the end of inflation also deflnes the surface 
on which inflation ends, which in this case is an ellipse. Varying /, g and h determines the 
orientation and the eccentricity of the ellipse in the {</), x} space; with g = deflning an ellipse 
aligned with the axes, and h = f with g = deflning a circle. With this in mind, to determine 
which contribution in Eq. (6.3) dominates we look at the relationship between /, h and g. 
For g ^ f,h the Ce term dominates over Cinf [60, 61]. We consider this case in the following 
subsection, and in the remainder of this section focus on the scenario where C — Cinf where 
the spectrum, bispectrum and trispectrum will evolve as described in our previous analysis 
for a two-field model, and there is no significant additional contribution to C at the end of 
infiation. The characteristics and various observational predictions of this model have been 
studied extensively in the literature, see for example Refs. [26, 27, 38, 39, 49, 52, 57, 62-69]. 

We restrict our attention to the case in which r/^ and rjy, are both positive, and so the 
potential during infiation has a valley-like shape. This provides an illustrative example of 
our discussion in §5.2. If either r]^ or rj^ were negative, the potential would contain a ridge. 
We further assume that ijfj, > rj^, and that the initial conditions are such that < x* *^ Mpi, 
whilst (p* ~ 0{Mpi) is significantly displaced from zero. Hence, the trajectory initially rolls 
towards the minimum along the 4> axis before turning and rolling very slowly towards the x 
axis. 

These initial conditions imply a small 9* , and since the potential is vacuum-dominated, 
a c:i 9 and the heatmaps are particularly simple to interpret. The initial (horizon crossing) 
condition corresponds to a position on the diagonal in the bottom left hand corner of any 
given heatmap. As the trajectory begins its slow turn into the valley, the point on the 
heat-map proceeds vertically upwards and we can see that a strong non-Gaussian signal is 
expected during the early stages of the turn. When the angle grows larger, the trajectory 
moves out of the regions of the heatmaps where / and r are large and so the non-Gaussianity 
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Figure 6. Evolution of sign(/NL)(6/5/NL)^ (solid blue line), tnl (dashed green line), and ^nl (dot- 
dashed red line) for the two-field hybrid model described in the text, calculated using the analytic 
formulae presented in the appendix. The vertical dashed line at 60 efolds indicates where we choose 
the waterfall field to terminate inflation. We also continue the evolution past this point, as if the 
waterfall field were absent, to show the subsequent peak and decay of the non-Gaussianity. 

decays. Since the large non-Gaussianity occurs during the early stages of the turn we see 
that 7]ss ~ % is a constant whilst all components of the .^?-^ parameter are zero. 

In Fig. 6 we show the evolution of /nli "^nl and gfNL for the potential (6.1). The 
parameter values used were r/^ = 0.09 and r]^ = 0.0025 with the initial conditions taken as 
(p* = 0.9Mpi and x* = O.OOlMpi. The evolution of the trispectrum is exactly as we expect: 
/nl and tnl both begin at negligible values, grow to large positive peaks and then decay again 
to negligible values. We confirm that the peak amplitudes are as expected from Eq. (5.4) in 
§5.2 to within 10%. Furthermore, we see that tnl peaks before /nl in agreement with the 
discussion regarding Eq. (3.16) and also in agreement with the results in §5.2. This quadratic 
shape of the potential gives us a small q'nl as discussed in §4.2.2. 

After 60 efolds, where we assume the evolution is terminated by the waterfall field, we 
find rig = 1.07, which is clearly in violation of the bounds on the spectral index placed by 
WMAP, which demands a red spectral tilt for priors of zero running and tensor to scalar 
ratio. As this model predicts neither a significant tensor to scalar ratio (> 10^^) nor a large 
negative running, we do not consider this model a viable candidate for inflation. Moreover, 
we know that there does not exist a choice of parameters or initial conditions (where both 
fields are less than the Planck scale) which can support a red tilt [26]. Nevertheless, we still 
consider the model a useful illustration of our analysis, and in particular of the a potential 
with a valley feature. It is also of note that a full treatment of the hybrid model may lead to 
the production of cosmic strings (for example see Refs. [70, 71]) and if cosmic strings were 
to contribute fractionally to the temperature power spectrum, then the best fit parameters 
would support a blue spectral tilt (see for example Refs. [72, 73]). 

6.1.1 The role of the end of inflation hypersurface 

As we have described, in addition to allowing for interesting phase space evolution during 
inflation, multi-field models of inflation can also have an end of inflation surface different 
from one of uniform energy density [57, 60, 61, 68, 74, 75], and this leads to the the extra 
contribution to the curvature perturbation ^e discussed above. If this extra contribution 
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dominates then the overaU curvature perturbation can be significantly non-Gaussian and so 
the trispectrum may be enhanced. 

Here we do not attempt a full analysis of the statistics when ^e is accounted for, but 
rather offer one simple example. We assume Ce ^ Cinf; ^-nd consider the case where x is held 
at zero during inflation and the parameters r]^ and r]^ are equal. The slow-roll parameters 
in this case are given as [60] 

2 

and fl'NL and tnl parameters at the end of inflation are found to be 

50(7 
5NL - -r— (3/nl iVss - 2e) + r/ss) , 
54 

/NL = % f ^ - :^V (6.5) 



Sg \9 2/ 

where all quantities are evaluated at the end of inflation. 

We see that tnl closely follows the single-field result tnl — (6/nl/5)^ and so may 
be detectable if /nl is observably large as found in ref. [60]. When /nl ^ 1, we see that 
9nl oc S'/nl X C'(e) and when /nl ^ 1 then (7nl oc g x 0{e). In both cases ^nl is negligible 
unless one considers a very large coupling g that is inconsistent with effective field theory 
and so we see that, for this particular model, modifying the end of inflation hypersurface 
produces non-Gaussianitites that are accounted for purely by the bispectrum. 

In Ref. [64], the author relates the non-Gaussianity predictions from the end of inflation 
hypersurface to the geometry of this surface. The author finds ^nl niay be significant when 
the curvature changes rapidly along the surface. Our findings here agree with this conclusion, 
because the parameter choices {g ^ 1) which lead to an enhanced gfNL correspond to a large 
ellipticity and have the trajectory lying very close to the semi-major axis. 

6.2 Axion Potential 

The model with the potential 

W=^g4>' + A'{l-cos{27rx/f)), (6.6) 

is inspired by the N-axion model studied in Refs. [51, 53]. It is similar to the model introduced 
in Ref. [38] as an example of a two-field potential which has an interesting transitory evolution 
as well as a large asymptotic value of non-Gaussianity at the adiabatic limit. Here we have 
chosen the (p field to self-interact through a quartic potential, which will allow a choice 
of parameters and initial conditions for which C can become constant during the slow-roll 
evolution, while still giving a value of the spectral index consistent with observations^. The 
potential (6.6) contains an infinite number of successive ridges and valleys. Taking initial 



^In Ref. [39], an example was given in which a large non-Gaussianity was reached during slow-roll inflation. 
We note, however, that the spectral index of this example was much too low to be consistent with observation. 
By choosing a quartic potential in this work, our aim is to show that there is no barrier in principal to having 
a model which reaches the adiabatic limit during slow-roll inflation with a consistent value of Us . 
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conditions such that x* is sufficiently close to //2, while (j) starts its evolution at a large 
field value {(jf ~ 23Mpi is required for roughly 60 efolds), the trajectory is initially almost 
entirely in the (j) direction and rolls along the ridge defined by x = //2- The initial angle 9* 
is therefore close to zero, and as the trajectory slowly turns off the ridge, 9 increases. 

From the the heatmaps, one can see that the initial growth of the angle 9 will lead to 
an increase in the magnitude of /nl and tnl as a passes up through the 'hot-spot' region 
where /nl and tnl peak. Since the third derivative of the potential is insignificant in the 
early stages of this turn, ^^nl is not significant. After the initial turn from the ridge-top, 
the trajectory will fall quickly down the steep side of the ridge, before turning back towards 
its original direction as it reaches the valley bottom. This leads to a decrease of a back 
through the hot-spot region once more. Importantly, the potential W drops by a small (but 
non- vanishing) amount between these two turns. From Eq.(4.2) this means that the initial 
growth of a is slightly larger than the subsequent decay and so a will not make it to zero 
as ^ — )• 0. Instead, it will end up with a small positive value. We therefore expect the non- 
Gaussianity parameters to take large constant values asymptotically, where we also know 
that the phase space trajectory ends up in a valley and so the asymptotic formulae for q'nl 
in Eq. (4.5) will be an accurate prediction. 

In Fig. 7 we give the evolution of /nl, ''"nl and (^nl for a specific realisation of this case 
described above. The parameter values are / = IMpi and h.^/g = (25/27r)^Mp[, with the 
overall normalisation fixed to agree with WMAP7 power spectrum amplitude. The initial 
conditions are (jf = 22.5Mpi and x* = //2 — O.OOlMpi which gives us the HCA value of a as 
0.022 which is small in agreement with the above discussion. In this case, the final constant 
spectral index has a value of ris = 0.949 which is within the WMAP7 95% contours. The 
evolution is exactly as we expect, with /nl and tnl both beginning with negligible values 
and /nl then evolving to a large negative peak while tnl grows to a large positive one, 
before both peaks decay. Despite the potential dropping quite significantly before the axion 
rolls, the peak values of /nl and tnl are described by the formulae in §5.1 with about 30% 
accuracy. The subsequent evolution is also interesting. As the trajectory evolves into the 
valley there is a positive spike in /nl and (7nl typical of this evolution. Furthermore, g'NL 
also increases significantly as we move away from the Vt = 1 regime and we see that there 
is a delay between the evolution of /nl and tnl as discussed in §5.1. Finally, the constant 
asymptotic values of /nl, ''"nl and /^nl are reached where we find Eq. (4.5) holds very well. 
In this case tnl ~ (6/5/nl)^) though this is not a necessary consequence of reaching the 
asymptotic regime, unlike the relation between tnl and (7nl- 

6.3 Inflection point model 

The discussion in §4.2.2 motivated sloping inflection features in the potential as a scenario 
in which the condition (7nl S> tnl might be found in an adiabatic regime where Q — ;■ 0. We 
now present such a model, but note that the tuning required to find parameters such that 
5fNL is both large and dominant in the adiabatic limit is severe. We expect, therefore, that 
the simple relation (4.5) will hold for the vast majority of models. 

We choose an infiection potential similar to that studied in Ref. [39] which takes the 
form 

W = Wo (^</.^ + Vo + hx+ ^Ax^ + ^MX' ) . (6.7) 

The value of Wq is fixed by the WMAP power spectrum, and there is an inflection point 
feature at x = 0. We further assume that Vq and // are fixed by the requirement that there 
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Figure 7. Evolution of sign(/NL)(6/5/NL)^ (solid blue line), tnl (dashed green line), and ^nl (dot- 
dashed red line) for the axion model discussed in the text, calculated using the analytic formulae 
presented in the appendix. 

is a minimum of the potential at Xmin = —p, where p is taken to be positive. We note that 
beginning some way above/below the inflection point leads to a negative/positive asymptotic 
value of /nl as explored in Ref. [39] . 

From an exploration of the parameter space for this model we have found it very difficult 
to obtain a large and dominant ^'nl whilst satisfying the constraints of slow-roll. We present 
a compromise scenario where the x field rolls when e ~ 0.3. Such a model has parameter 
values h = O.OSM^j, A = lO'^Mpi and p = O.lMpi and the evolution begins with (/>* = 22.5Mpi 
and X* = ^- We compare the evolution from our analytic expressions to the numerical 
non-slow-roll evolution calculated using a finite-difference code identical to that used in Ref 
[38]. In the non-slow- roll case, the scalar field velocity perturbations decouple from the field 
perturbations and so the phase space of the dynamics is enlarged. We see the effect of this by 
the presence of oscillatory behaviour in the non-Gaussianity parameters and in the spectral 
index. When the deviation from slow-roll is not too large, the non-slow-roll evolution mimics 
a superposition of the the slow-roll evolution and these rapid oscillations. We therefore find 
that our analytic results are surprisingly applicable when slow-roll begins to break. For this 
example, the analytic spectral index is Ug = 0.964 and this is within 1% of the non-slow-roll 
value. Our calculation gives (7nl = —432 which agrees with the non-slow-roll code to within a 
factor of two. In general, if the evolution of the observables for a given model have not settled 
down to a constant value before the end of slow-roll inflation then one must use numerical 
methods to obtain an accurate result. However, a crude estimate is very easy to calculate 
through the slow-roll analysis. To conclude, in this final example we have tried to engineer 
a model to break from the simple 5'nl~''"nl result in Eq. (4.5) but we have found that this is 
very hard to achieve. 

® A realistic application of our formulae requires that the model is insensitive to changes in the horizon 
crossing conditions within the range of the quantum scatter. Such changes do not affect the predictions of 
this model. 
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7 Conclusions 

In this work we have studied non-Gaussianities in two-field slow-roh inflation. By extending 
and simpUfying the 'heatniap' analysis of Byrnes et al. [27] we have shown, for separable po- 
tentials, the regions of the parameter space that can give rise to large non-Gaussianities, both 
in the bi- and tri-spectra. This analysis also benefits from showing the explicit dependence of 
these non-Gaussianites on the shape of the inflationary potential, and the Horizon Crossing 
Approximation can be easily applied to our results. The analyses in the paper allows us to 
make the following conclusions: 

• We have found that the regions of the heatmap parameter space that produce large 
values of the bispectrum parameter /nl are also capable of producing large values of the 
trispectrum parameters tnl or s^nl- Our results confirm that a necessary requirement 
for a large local non-Gaussianity is that the horizon crossing field velocities must be 
dominated by one of the two fields. 

• The relationship between tnl and /nl is shown explicitly for separable potentials and 
is given in Eq. (3.16). 

• We have found that in the adiabatic limit the non-Gaussianity parameters gfNL and tnl 
can be related by Eq. (4.5), which is confirmed for a variety of models considered here. 
This relation can be modified when the isocurvature third slow-roll parameter ^1^^ is 
significant. 

• We have shown that the ridge and valley features in the inflationary potential can lead 
to large tnl and /nl parameters, with tnl being the flrst to peak. g'NL can be large if 
the potential is described by terms beyond quadratic order. It is possible, but very hard 
to engineer, that g'NL is the dominant statistic if inflation reaches a natural adiabatic 
limit. 

• We have employed the heatmaps in a dynamic way, and shown that they can be used as 
qualitative tools, in order to understand and predict the evolution of non-Gaussianities 
in different models of inflation. 

Forthcoming Planck data will further constrain primordial non-Gaussianities. Since our 
analysis relates these observations to the shape of the multi-field inflationary potential, it may 
be that if a signal is detected one could infer information about the inflationary potential 
from observation. This is an exciting possibility. Moreover, the non-Gaussianity parameters 
considered here are not the only parameters that can act as signatures for specific shapes 
in the potential. In particular, it would be valuable to extend our analysis to consider the 
running of the spectral index and the running of /nl? ffNL and tnl, which we intend to return 
to in future work. 
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A Summary of results for two-field SN calculations. 

A.l Sum-separable potentials 

This section follows from the work of Vernizzi and Wands [25] and considers a potential with 
the sum-separable form W = U{(/)) + V{x)- The first derivatives of A^ are 

U* + Z 



MpiN^^ = ■ 


^jn' 


u - 


^Pi^.x = ■ 


V 




V^^x' 
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._Ve^- 


-Ue^ 



w* 
V* - z 



(A.l) 



w* 

V A 

where 

(A.2) 
e 

For notational ease we have dropped the 'c' label that is commonly attached to quantities 

evaluated on a later-time uniform density hypersurface. For the purposes of calculating 

the trispectrum we shall also need to decompose rnj and ^?-^, into the kinematic basis on the 

uniform density hypersurface. In the original {(p, x} frame the potential is sum-separable and 

the only non-zero values of r/jj and ^J-^ are those for which all of the indices are identical, 

which allows us to use the single-index notation 77^ and A. After rotating into the kinematic 

basis with fields a and s, the potential generally loses it separable form and so it is necessary 

to use all of the indices. The three r/ components are 

'/(TCT ) 



\fWxt ^ 

^^s = — ^ — {r]^-r]4,), 

rjss = ^i±^^ . (A.3) 



The components of ^"^.^^ in the kinematic frame are 



,3/2 c2 _ 3/2 t2 , 3/2 ^2 

3/2 C2 = ^_ ATT C2 _ /— /:2 



.3/2^2 =,,_/e^^2^ /^^2 



^ass ~ '^0V ^X ^x ' X 

e"'esss = eTil-ellHl- (A.4) 

The second derivatives of A^ follow by differentiation of Eqs. (A.l), giving 
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where 



W* ]j 2 d(j)* W* \l 2 dx* 

Taking the next derivative we find 
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The other derivative of A follows by a simple permutation of Eq. (A.12) under the joint 
permutations {u v} and {</> -H- x} which has the effect of negating B^ whilst A remains 
unchanged. We take these results together to find 
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We now rewrite the horizon crossing slow-roll parameters in terms of their kinematic 
counterparts such as r/*g. There are three kinematic rj parameters and four ^^ parameters 
which means that there is no unique way in which to do this. We shall choose to rewrite r/^ 
and r/* in terms of r]*^ and 77*^ and ^1 and ^* in terms of Qg^ and S,ass • This gives the 
relations 



V4> = Vs 



Vx ~ ^ss 
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Substituting Eq. (A.15) into Eq. (A. 14) and simplifying we find 
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A. 2 Product-separable potentials 



For a potential with the two-field product-separable form W = U{(f))V{x): the first derivatives 
of N are 
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In the {(f), x} frame the potential is product-separable which means that any of the parameters 
rjij or ^f-j^ with mixed derivatives can be written in terms of lower-order slow-roll parameters 
that do not have mixed derivatives. This again allows us to use the single-index notation r/^ 
and ^'i for the remaining terms. In the kinematic basis the three t] components are 

'ha — ) 
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The components of the ^f-f^ tensor are 
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e'/' ds = ef il - 4^ il + 6V^e^(ex^<^ " e<^^x) • (^-19) 

The second derivatives of A^ follow by differentiation of Eqs. (A.17), giving 
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where we have substituted 



Ap = uv'qss-, (A.21) 
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to put these equations into a form similar to that found for the sum-separable potential. 
Taking the next derivative we find 
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As with the sum-separable case, Ap is symmetric under switching the fields whilst Bp is 
anti-symmetric. We take these results together to find 
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One can rewrite the horizon crossing slow-roll parameters in terms of their kinematic 
counterparts as 
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Substituting Eq. (A. 29) into Eq. (A. 28) and simplifying we find 
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B Simplification of Trispectrum expressions 

The function \tq\, as defined in Eq. (3.13), is bounded well within our limit of ten and so 
this term may be immediately neglected. Furthermore, we can use Eq. (3.6) to manipulate 
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various other terms. In the sum-separable case we find 
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and we find that all of the functions in square brackets never have a magnitude greater 
than ten and so these terms represent variations in the trispectrum that are significantly 
smaller than observables will ever probe. One can easily check that these results follow 
analogously for product-separable potentials and so these terms may be neglected for both 
types of separable potential. 

Analogously to Eq. (3.6) there exist formulae relating the various ^?-^ components. For 
sum-separable potentials we have 
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whereas for product-separable potentials these relations are of the form 
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Using Eqs. (B.2) and (B.3) we can then simplify the full expression for ^nl in Eq. (3.10) as 
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where we have defined 54 = ^ (rs sin 29* cos 29* — T2) in the expressions above. The term in 
square brackets in the last line can never be large and so is neglected. 

The product-separable case follows similarly, yielding the same answer, however the 
calculation is unsurprisingly more involved. Manipulating three of the terms in Eq. (3.12) 
for (7nl by using Eq. (B.4) we find 
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We now expand the last term of Eq. (B.7) by substituting for ^*g.^ using Eq. (B.5). We also 
use Eq. (3.6) to rewrite r]*^ as sin20*7/*^ + (1 — sm29*)^ ta.n29* {rjss — "q^u + 2e) and so we 
ultimately find 






-T3sin2 26l* 
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-e*2 + 3e*C cos2r + ^(1 - sin2r)(C - r?;, + 2e*)j . (B.8) 

The term in square brackets is always negligible and this multiplies terms no larger than 
0{e* ) and so may be ignored, leaving the same simple result that we found for sum-separable 
potentials. 

After these various terms have been neglected from Eqs. (3.9) to (3.12), we can simplify 
the remaining terms by rewriting them by means of the trigonometric relations 

n = T + 2/ + 1 , 

r4 = 2/i(l + /), 
Tj = 4(r + /) , 
2g2 = T-f. (B.9) 
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